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COMPARISON  OF  DETECTORS  IN  THE  PRESENCE  OF  SIDELOBE  JAMMING 


INTRODUCTION 

Automatic  detection  systems  have  been  designed  with  the  philosophy  that  they  will 
suffer  a  small  detection-sensitivity  loss  in  thermal  noise  and  at  the  same  time  the  number  of 
undesirable  detections  (i.e.,  noise  false  alarms  and  clutter  residues)  will  be  limited  so  that  the 
automatic  tracking  system  will  not  be  overloaded.  In  this  report,  the  design  of  an  automatic 
detector  will  be  considered  for  the  case  of  sidelobe  noise  jamming.  The  basic  idea  is  to  con¬ 
figure  the  detector  so  that  targets  can  be  detected  when  the  jammer  is  in  sidelobe  nulls  of 
the  antenna. 


MATHEMATICAL  FORMULATION 

If  the  jamming  is  white  Gaussian  noise,  the  density  of  the  ith  output  pulse  x,  from  an 
envelope  detector  is 


PiXi^i)  =  +A?)t2°f\  *o 


(1) 


where  A,-  is  the  signal  amplitude  and  the  signal-to-noiae  ratio  (S/N)  is  10  log  (Af/2of).  The 
optimal  detector  is  the  likelihood  ratio 


n 

L  -  *  0)]  ,  (2) 


where  , ....  xn  are  n  independent  samples.  A  target  is  declared  when  L  is  greater  than  the 
threshold  which  determines  the  probability  of  false  alarm  (Pfa).  Using  a  small  signal  approxi¬ 
mation,  we  can  show  that  the  likelihood  test  is  equivalent  to  comparing 


(3) 
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to  an  appropriate  threshold.  If  the  noise  samples  are  identically  distributed  (oy  =  o  for  all  i) 
and  the  {A,}  are  set  to  one,  then  (3)  reduces  to 


X>.! 


?  >Ko2 


Since  a2  is  usually  unknown,  a2  must  be  estimated.  The  implementation  of  (4)  is  the  cell¬ 
averaging  CFAR  (constant  false  alarm  rate)  detector  discussed  by  Finn  and  Johnson  [1J  and 
shown  in  Fig.  1.  This  detector  uses  the  shift  registers  (SR)  to  save  the  last  N  samples  in  each 
range  cell,  integrates  the  samples  in  a  moving  window  (MW),  and  then  performs  the 
thresholding  (CFAR)  function. 

To  investigate  the  effects  of  jamming  on  this  detector,  let  us  consider  the  return  signal 
and  jammer  powers,  Pg  and  Pj  respectively.  If  we  define  S  by 

Ps  =  SG*  ,  (5) 


where  Gs  is  the  radar-antenna  power  gain  (one  way)  in  the  direction  of  the  target  divided  by 
the  midbeam  gain,  and  we  define  J  by 


-  Cell-average  CFAR  detector:  SR  is  a  ahift  register, 
MW  is  a  moving  window,  and  C  is  a  comparator 
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PJ=JGj 


(6) 


where  Gj  is  the  radar-antenna  power  gain  in  the  direction  of  the  jammer  divided  by  the  mid¬ 
beam  gain,  S  and  J  represent  received  signal  and  jamming  powers  independent  of  the  an¬ 
tenna  geometry  with  respect  to  the  target  and  jammer.  Using  these  normalized  powers  and 
letting  Pn  be  the  received  thermal  noise  power,  we  can  write 


S 

N 


SG2s 

Pn  +  JGj 


(7) 


Consider  a  hypothetical  200  n.  mi.  radar  that  has  16  pulses  on  target.  From  Robertson 
[2] ,  a  S/N  =  4  dB  per  pulse  is  required  to  achieve  a  probability  of  detection  (PD)  of  0.9  and 
a  pfa  of  10-6.  Thus,  at  a  range  of  100  n.  mi.,  S/N  =  16  dB  per  pulse,  a  value  large  enough  so 
that  detections  can  be  declared  on  a  single  pulse.  Now  consider  the  target-jammer  separation 
with  respect  to  the  sin  0/0  antenna  pattern  shown  in  Fig.  2.  The  jammer  is  in  a  null  of  the 
antenna  pattern,  and  the  signal  power  is  only  reduced  a  few  dB;  i.e.,  Gg>  Gj.  Conse¬ 
quently,  as  shown  in  Fig.  3,  the  S/N  (given  by  (7))  is  large  for  one  of  the  pulses  as  the 
antenna  beam  sweeps  by  the  target,  and  theoretically  the  target  can  be  detected.  However, 
if  the  cell-averaging  CFAR  detector  is  used,  the  target  will  not  be  detected  because  the 
average  signal  power  is  small  in  comparison  to  the  average  jamming  power. 

For  example,  if  the  geometry  in  Fig.  2  is  used,  if  S/Pn  =  20  dB,  and  if  S/J  =  -20  dB, 
the  S/N  ratios  in  Table  1  are  obtained.  Examination  of  Table  1  reveals  that  the  highest 
S/N  (13.6  dB)  is  obtained  when  the  target  is  located  approximately  8  pulses  before  the  beam 
center.  If  no  integration  were  performed  (i.e.,  if  the  decision  were  made  on  a  single-pulse 
basis),  a  S/N  =  13.6  dB  would  yield  a  PD  >  0.9.  However,  if  the  cell-average  CFAR  detector 
is  used,  the  equivalent  input  S/N  is  given  by 


(8) 


For  the  data  in  Table  1,  the  greatest  ( S/N)c  ,  is  0.38  (-4.2  dB),  and  this  value  is  obtained 
by  integrating  the  16  pulses  from  - 14  to  1.  From  Robertson  [2]  (assuming  there  is  no  dif¬ 
ference  between  linear  and  square-law  detectors),  the  PD  corresponding  to  n  «  16  and 
S/N  -  -  4.2  dB  is  less  than  0.05. 
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Fig.  3  —  Signal-to-noiie  ratio  per  pulae  aa  antenna  sweep*  over  target 
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Table  1  -  S/N  Per  Pulse 


Pulse  Position 
w.r.t.  Beam  Center 

S/N 

(dB) 

-15 

-21.5 

-14 

-17.0 

-13 

-14.0 

-12 

-10.1 

-11 

-  6.1 

-10 

-  1.2 

-  9 

4.4 

-  8 

13.6 

-  7 

6.4 

-  6 

2.8 

-  5 

0.9 

-  4 

-  1.1 

-  3 

-  2.5 

-  2 

-  4.0 

-  1 

-  4.5 

0 

-  5.0 

1 

-  6.0 

2 

-  7.0 

3 

-  8.0 

4 

-  9.0 

5 

-10.0 

On  the  other  hand,  if  the  test  cell  is  normalized  (divided)  by  the  average  of  the  reference 
cells  and  then  integrated  as  in  Fig.  4,  the  effective  S/N  for  this  ratio  detector  is  given  by 


(S\  _ ^'2 _ 

W*  *  nl,  Pn{i)  +  PjV) 


For  the  data  in  Table  1,  the  greatest  ( S/N )R  is  2.32  (3.7  dB);  this  value  is  obtained  by  in¬ 
tegrating  the  16  pulses  from  -11  to  4.  If  Robertson’s  curves  [2]  are  used  (although  they 
were  derived  for  a  linear  detector,  they  will  yield  approximate  results  for  the  ratio  detec¬ 
tor),  the  Pl}  =  0.9. 

The  reason  the  ratio  detector  is  very  effective  can  be  understood  if  we  reexamine  (3): 
since  o,  is  unknown  and  varying  from  pulse  to  pulse,  o,  cannot  be  placed  on  the  right-hand 
side  of  the  inequality  as  in  (4).  It  can  readily  be  shown  that  no  uniformly  most  powerful 
test  exists  (i.e.,  the  optimal  test  depends  on  specific  values  of  {o, )  ).  A  logical  way  of  pro- 
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Fig.  4  --  Ratio  detector:  SR  is  a  shift  register,  MW  is  a 
moving  window,  NR  is  the  number  of  reference  cells, 
and  C  is  a  comparator 


ceeding  is  to  use  the  generalized  likelihood  ratio  test  of  Van  Trees  [3] ,  replacing  the  un¬ 
known  parameters  with  an  estimate.  Thus,  it  can  be  seen  that  Fig.  4  is  just  an  implementa¬ 
tion  of  the  generalized  likelihood  ratio  detector.  The  detector  output  is 


E 


xfV) 


k*  1 


(10) 


where  £,(/)  is  the  tth  envelope  detected  pulse  in  the  /th  range  cell  and  NR  =  2m  is  the  num¬ 
ber  of  reference  cells.  The  denominator  is  an  estimate  of  of  and,  essentially,  the  detector 
sums  signal-to-noise  power  ratios. 

It  should  be  noted  that  the  square-law  detector  should  be  used  rather  than  the  linear 
detector,  since  the  problem  is  equivalent  to  the  collapsing  loss  problem.  However,  instead  of 
the  signal  disappearing  for  many  pulses,  the  noise  (jamming)  increases  for  many  pulses. 

Since  it  has  been  shown  by  Trunk  (4)  that  the  collapsing  loss  is  much  smaller  for  a  square- 
law  detector,  this  detector  should  be  used. 
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DETECTOR  CALCULATIONS 

To  compare  the  performance  of  the  ratio  detector  with  that  of  the  cell-averaging 
CFAR,  we  generated  PD  curves  using  a  Monte  Carlo  simulation.  First,  by  use  of  the  concept 
of  importance  sampling  (see  appendix),  threshold  curves  were  generated;  they  are  shown  in 
Fig.  5.  The  PD  curves  for  the  two  detectors  in  thermal  noise  were  found  by  a  straight¬ 
forward  simulation  and  are  shown  in  Figs.  6  and  7.  The  ratio  detector  requires  0.2  dB  more 
signal  to  yield  the  same  PD  as  the  cell-average  CFAR. 

The  performance  of  the  detectors  in  sidelobe  jamming  is  quite  complicated;  it  is  a  func¬ 
tion  of  S/P„ ,  S/J,  and  the  angular  separation  of  target  and  jammer.  Curves  of  PD  versus  tar¬ 
get-jammer  separation  in  beamwidths  for  S/J  =  -  20  dB  and  S/Pn  =  16  and  20  dB  are  shown  in 
Fig.  8.  We  calculate  the  PD  by  finding  the  maximum  detector  outputs  as  the  beam  sweeps 
over  the  target  and  then  comparing  these  maximum  values  to  the  appropriate  thresholds.  In 
all  cases,  the  performance  of  the  ratio  detector  is  better  than  that  of  the  cell-average  CFAR. 
Also,  the  greater  the  value  of  S/Pn ,  the  greater  the  performance  difference. 


Fig.  5  —  Threshold  values  for  cell-average  CFAR 
and  ratio  detectors 
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S/N=I6  dB/ pulse  S/N  =  20  dB/pulse 


TARGET- JAMMER  SEPARATION  (BEAMWIDTHS) 


Fig.  8  —  Probability -of-detection  curves  for  detectors  in  sidelobe  jamming: 
SIJ  =  -  20  dB  and  Pfa  =  10~6 


ANGULAR  ESTIMATES 

When  a  target  is  detected  with  the  ratio  detector  (with  the  jamming  being  in  an  antenna 
null),  it  is  impossible  to  obtain  an  accurate  azimuth  estimate.  Specifically,  for  very  high  S/Pn 
and  very  low  S/J  the  azimuth  error  is  uniformly  distributed  in  an  interval  approximating  the 
beamwidth.  Also,  the  azimuth  estimates  are  not  independent  from  scan  to  scan,  because  the 
estimate  is  essentially  the  midbeam  antenna  position  when  the  jammer  is  in  a  null  and  the 
target  is  in  the  main  beam.  Thus,  the  azimuth  estimates  follow  the  jammer  movement  (with 
a  large  bias)  until  the  target-jammer  separation  changes  so  that  the  target  is  detected  when 
the  jammer  is  in  a  different  antenna  null.  At  this  time,  the  estimate  will  change  by  a  beam- 
width. 

The  type  of  behavior  one  obtains  is  illustrated  in  Fig.  9.  In  this  example,  the  jammer  is 
initially  at  90°  and  is  moving  0.3° /scan,  the  target  is  initially  at  87°  and  is  moving 
-0.4° /scan,  and  the  nulls  are  separated  by  3°.  The  estimated  azimuth  lags  the  jammer  by 
one  beamwidth  initially,  later  by  two  beamwidths,  and  finally  by  three  beamwidths.  Special 
care  must  be  taken  to  incorporate  this  type  of  behavior  into  the  tracking  algorithm,  since 
the  angle  estimates  are  uniformly  distributed  about  the  true  azimuth  and  will  occasionally 
have  discontinuities  of  the  order  of  a  beamwidth,  whereas  normally  the  azimuth  estimates 
are  assumed  to  be  Gaussian  distributed. 


EXPERIMENTAL  RESULTS 

Theoretically,  the  performance  of  the  ratio  detector  is  much  better  than  that  of  the 
cell-average  CFAR.  Consequently,  it  was  decided  to  compare  these  detectors  using  recorded 
data.  Specifically,  two  S-band  TWTs  (located  on  building  #5  and  on  the  north  gatehouse 
at  bearings  of  approximately  167°  and  185°,  respectively,  from  the  radar)  were  used  to  jam 
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the  SPS-39  radar  at  CBD.  A  recording  system  recorded  65532  eight-bit  numbers  corre¬ 
sponding  to  a  designated  R-6  region  as  the  radar  rotated.  After  the  65532  data  words  were 
accumulated,  the  data  were  dumped  on  to  magnetic  tape  in  approximately  1.5  s.  Thus,  the 
system  was  ready  to  record  the  next  scan  of  the  designated  R-6  region  by  the  time  the  beam 
next  entered  the  region. 

The  data  presented  in  this  paper  were  taken  on  May  28, 1980.  Data  were  recorded 
from  two  different  R-6  regions.  The  first  region  started  at  a  range  of  50  n.  mi.  and  an 
azimuth  of  115°  and  consisted  of  64  range  cells  for  each  of  1024  transmitted  pulses.  The 
second  region  again  started  at  a  range  of  50  n.  mi.  and  an  azimuth  of  115°  but  consisted  of 
128  range  cells  for  each  of  512  transmitted  pulses.  Since  the  radar  PRF  was  500  pps  and  the 
radar  scan  rate  was  8  s,  the  first  region  covered  90°  of  azimuth  and  the  second  covered  45°. 
For  each  region,  three  different  sets  of  jamming  data  were  recorded: 

•  The  jammer  at  the  north  gatehouse  was  set  at  maximum  attenuation  and  the  enve¬ 
lope-detected  samples  were  recorded  as  the  jamming  power  (on  building  #5)  was  reduced 
in  6-dB  steps. 

•  The  jammer  at  building  #5  was  set  at  maximum  attenuation  and  data  were  recorded 
as  the  jamming  power  at  the  north  gatehouse  was  reduced  in  6-dB  steps. 

•  The  jammers  were  set  to  approximately  the  same  power  and  data  were  recorded  as 
both  jammers  were  attenuated  in  6-dB  increments. 
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Some  typical  recorded  data  from  region  1  are  plotted  in  Fig.  10  and  Figs.  11a  to  11/, 
where  the  returned  power  is  found  by  averaging  the  64  range  samples  available.  Since  the 
radar  was  not  transmitting,  these  plots  represent  the  sum  of  the  average  thermal  noise  and 
jamming  powers  as  a  function  of  scanning  angle.  If,  in  the  denominator  of  Eq.  (7)  (i.e., 

Pn  +  JG j),  JGj  >  Pn,  the  plots  represent  the  antenna  gain  G j.  Thus,  by  noting  that  the 
average  noise  level  in  Fig.  10  is  18  dB,  we  can  infer  that  the  mainbeam  jamming  (with  no 
attenuation)  is  70  dB  above  the  thermal  noise.  Since  there  are  several  local  minima  15  to 
20  dB  below  adjacent  maxima,  one  would  expect  the  ratio  detector  to  have  a  higher  proba¬ 
bility  of  detection  than  the  cell-average  CFAR. 


To  compare  the  two  detectors,  signal  is  added  to  the  recorded  data  to  simulate  the 
presence  of  a  target.  For  example,  since  there  are  approximately  21  pulses  between  the  3-dB 
antenna  points,  signal  can  be  introduced  into  any  range  cell  by 


rf-  [~4  sin  K(j  - 

L  *0-0 


+  Zj  cos  6i 


+  (Zj  sin  0.)a 


where  A  is  the  signal  amplitude,  K  =  1.3916/10.5,  j  is  the  indexed  azimuth  location  of  the 
target,  and  zi  is  the  ith  recorded  pulse  (PRF)  in  the  range  cell.  The  signal  is  first  introduced 
into  the  ninth  range  cell.  Then,  this  signal  and  the  jamming  signals  in  the  eight  range  cells 
on  either  side  of  the  ninth  cell  are  put  into  the  cell-average  CFAR  (Fig.  1)  and  the  ratio 
detector  (Fig.  4)  and  any  detections  are  noted.  This  calculation  is  repeated  for  the  target  in 
the  13, 17,  21, ...,  49,  and  53  cells,  yielding  12  detection  opportunities.  The  signal  level 
A  for  each  azimuth  is  adaptively  adjusted  so  that  the  number  of  detections  made  by  the 
ratio  detector  is  between  9  and  11  for  region  l.*  Then  the  number  of  detections  made  by 
both  detectors  is  recorded;  this  is  plotted  in  Fig.  12.  A  comparison  of  the  two  detection 

•If  the  signal  strength  is  greater  than  255,  this  value  exceeds  the  dynamic  range  of  the  recorded  data  and 
consequently  the  number  of  detections  is  arbitrarily  set  to  zero  for  both  detectors. 


Fig.  10  —  Returned  noise  power  vs  radar  azimuth  angle 
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Fig.  11  —  Returned  jamming  plu<  noiae  power  va  radar  azimuth  angle.  The  jamming  attenuation*  for  building 
*5  and  the  north  gatehouae  (NGH)  are  given  for  each  case. 


r*.  11  --  Returned  jamming  plus  noise  power  n  radar  azimuth  angle  (Continued) 


(c)  Attenuation*  Bid*  #8  »  27  dB  and  NGH  -  12  dB 
1 2  -  Companaon  of  ratio  detector  and  otl-arcr **»  CFAR 
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curve*  shows  that  there  are  several  instances  where  the  probability  of  detection  of  the  ratio 
detector  is  over  80%,  whereas  that  of  the  cell-average  C FAR  is  zero.  Furthermore,  these 
large  differences  occur  when  there  is  a  very  sharp  and  deep  null  in  the  antenna  pattern. 

Table  2  was  generated  to  try  to  give  a  number  which  describes  the  advantage  of  the 
ratio  detector  over  the  cell-average  CFAR.  Here,  the  overall  probability  of  detection  arrived 
at  by  averaging  over  all  azimuths  indicates  less  than  a  1-dB  advantage  for  the  ratio  detector. 
However,  this  does  not  tell  the  whole  story.  If  one  refers  to  Fig.  13,  which  gives  the  histo¬ 
gram  of  the  number  of  detections  for  the  first  entry  in  Table  2,  one  sees  that  14%  of  the 
time  the  cell-average  CFAR  did  not  detect  the  target  even  once,  whereas  the  average  number 
of  detections  was  10  for  the  ratio  detector.  Huts,  14%  of  the  time  the  ratio  detector  was 
more  than  4  dB  better  than  the  cell-average  CFAR. 


Table  2 


Jamming  Attenuation  (dB) 

Detection  Probabilities 

Difference 

(dB) 

Building  #5 

North  Gatehouse 

Ratio 

C.A.  CFAR 

21 

Max 

0.834 

0.556 

Max 

6 

0.836 

0.642 

27 

12 

0.841 

0.686 

IKS 

Max 

Max 

0.797 

0.760 

Fig  13  —  Histogram  of  number  of  detection*  at  each 
azimuth  for  the  ratio  detector  and  cell-average  CFAR 


NRL  REPORT  8449 


CONCLUSIONS 

The  problem  of  designing  an  automatic  detector  in  a  sidelobe-jamming  environment 
has  been  considered,  and  it  has  been  shown  that  a  ratio  detector  has  better  performance 
than  the  cell-average  CFAR  detector.  The  performance  of  the  detectors  is  essentially  equiv¬ 
alent  in  thermal  noise.  However,  in  sidelobe  jamming  with  S/P  =  20  dB  and  S/J  =  -  20  dB, 
theoretically  the  Pd  for  the  ratio  detector  is  greater  than  0.99,  whereas  the  Pd  for  the  cell- 
average  CFAR  detector  ranges  between  0.02  and  0.42,  depending  on  target-jammer  angular 
separation.  Furthermore,  experiments  using  recorded  jamming  data  confirm  the  advantages 
of  the  ratio  detector.  Specifically,  although  the  advantage  of  the  ratio  detector  (averaged 
over  all  azimuths)  is  less  than  1  dB,  more  than  a  4-dB  advantage  occurs  about  14%  of  the 
time.  Also,  it  is  worth  noting  that  the  ratio  detector  will  have  a  performance  advantage  in 
jamming  when  pulse-to-pulse  frequency  diversity  is  employed.  This  is  because,  when  fre¬ 
quency  diversity  is  used,  the  antenna  gain  in  the  direction  of  the  jammer  will  change  and, 
correspondingly,  the  received  jamming  power  will  vary  from  pulse  to  pulse. 

The  main  disadvantage  of  the  ratio  detector  is  that  it  will  detect  single-pulse  inter¬ 
ference.  However,  one  can  probably  obtain  performance  equivalent  to  that  of  the  ratio 
detector  by  using  two  types  of  detectors,  the  normal  cell-averaging  CFAR  and  an  Af  out  of 
N  detector  where  M  -  N  =  2  regardless  of  the  number  of  pulses  in  the  beamwidth.  This 
possibility  will  be  investigated  in  future  work. 
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Appendix 

IMPORTANCE  SAMPLING 

Although  Monte  Carlo  simulations  have  been  used  for  many  years  to  calculate  PD 
curves,  they  have  not  been  used  to  calculate  Pfa  curves  because  of  the  enormous  number 
of  repetitions  usually  required,  approximately  10 IPfa .  However  this  difficulty  can  be  over¬ 
come  by  the  use  of  importance  sampling.  The  fundamental  principle  of  the  importance¬ 
sampling  technique  is  a  modification  of  the  probabilities  that  govern  the  outcome  of  the 
basic  experiment  of  the  simulation  so  that  the  event  of  interest  (the  false  alarm)  occurs 
more  frequently.  We  then  compensate  for  this  distortion  by  weighting  each  event  by  the 
ratio  of  the  probability  that  this  specific  event  would  have  occurred  if  the  true  probabil¬ 
ities  had  been  used  in  the  simulation  to  the  probability  that  this  same  event  would  occur 
with  the  distorted  probabilities.  Consequently,  by  proper  choice  of  the  distorted  proba¬ 
bilities  the  number  of  repetitions  can  be  reduced  greatly.  For  instance,  the  mean  of  a  func 
tion  {?(*)  is  given  by 


E{Q(x))  *  fq(x)dP{x), 


where  P(x)  is  the  distribution  of  x.  We  can  estimate  the  mean  of  Q(x)  by  selecting  M  inde¬ 
pendent  samples  x,  from  P(x)  and  associating  the  probability  1/Af  with  each  event.  The 
E  Q{x)  can  be  estimated  by 


M 


i-l 


(Al) 


The  importance-sampling  technique  uses  the  Radon-Nikodym  derivative  to  express  the 
mean  value  of  Q(x)  by 


F{0<x)}  -Jq(x) 


dP(x) 

dG(x) 


dC’(x) , 


where  (»'(* )  is  a  distribution  function.  We  can  now  estimate  the  mean  £{©(x)}  by  selecting 
Af  independent  samples  from  G(x)  and  associating  the  probability  dP[x,)lMdG(x,)  with  each 
event  ©<x, ).  Thus  E  {©(x ) J  is  estimated  by 
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M 


M 


E  Qtx'] 


i- 1 


dPjXj) 

dG(Xj) 


(A2) 


Since  (Al)  and  (A2)  are  both  unbiased  estimates  of  Q(x),  it  is  possible  to  select  G(x)  so  that 
the  variance  of  (A2)  is  less  than  the  variance  of  (Al). 

In  our  problem  of  determining  the  threshold  T  for  a  desired  Pfo  for  the  ratio  detector, 
it  is  necessary  to  estimate  the  distribution  curve 


P(Z<T)=1-  Pfa  ,  (A3) 


where 


n  xf(J) 

Z  =  E  — - 1 - ; 

imi  ^  jr  \xf(j + 1 + *) + xfo - 1  -  k) 

k~i 


(A4) 


and  the  x,(k)  are  independent  Rayleigh  variables.  The  straightforward  way  of  generating  the 
squared  samples  is  by 


xf(k)  =  o2(-  2  In  Uik )  for  all  k  , 


(A5) 


where  Uik  are  independent  random  numbers  uniformly  distributed  on  the  interval  (0,  1).  To 
estimate  (A3),  we  form  M  independent  ratios {Zj,  «/=  1, ...,  fusing  (A4)  and  (A5);  the  esti¬ 
mated  distribution  is 


P(Z  >  T) 


J- 1 


where 


1,  Zj>T, 


=  0,  Zj  <  T  . 
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Importance  sampling  differs  from  the  previous  procedure  by  generating  samples  using 


xf(j)  =  a 2(- 2  in  Uy)  (A6) 


and 


xf(k)~oH-  2\nVik)  k*j,  (A7) 

where  a  >  o,  a  technique  which  yields  more  false  alarms.  From  ( A4),  ( A6),  and  (A7),  : 

M  ratios  Z j  are  generated.  Then  the  estimated  distribution  is  j 


J*  1 


where 


ij»l,  Zj>Ty 


“0,  Zj<T , 


and 


